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Abstract 

We model the dynamics of asset prices and associated derivatives by consideration 
of the dynamics of the conditional probability density process for the value of an asset 
at some specified time in the future. In the case where the price process is driven by 
Brownian motion, an associated "master equation" for the dynamics of the conditional 
probability density is derived and expressed in integral form. By a "model" for the 
conditional density process we mean a solution to the master equation along with the 
specification of (a) the initial density, and (b) the volatility structure of the density. 
The volatility structure is assumed at any time and for each value of the argument of 
the density to be a functional of the history of the density up to that time. In practice 
one specifies the functional modulo sufficient parametric freedom to allow for the input 
of additional option data apart from that implicit in the initial density. The scheme 
is sufficiently flexible to allow for the input of various types of data depending on the 
nature of the options market and the class of valuation problem being undertaken. 
Various examples are studied in detail, with exact solutions provided in some cases. 
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1 Introduction 

This paper is concerned with modelling the dynamics of the volatility surface. The problem 
is of great practical interest to traders, and as a consequence has an extensive mathematical 
literature associated with it. In this brief report, we shall not attempt to review earlier 
work in the area, but refer the reader, for example, to Schonbucher (1999), Gatheral (2006), 
Schweizer & Wissel (2008a,b), Carmona & Nadtochiy (2009, 2011), and references therein. 
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Put informally, the general idea of the paper is as follows. We fix a numeraire, and write {A t } 
for the value process of some tradable financial asset expressed in units of that numeraire. 
We fix a time T and assume that no dividends are paid from time up to T. Letting Q 
denote the martingale measure associated with the chosen numeraire, we have 

A t = E®[A T \F t ], (1.1) 

for t <T, and 

C tT (K)=E Q [(A T -K) + \F t ], (1.2) 

where Ctr{K) denotes the price at time t of a T- maturity call option with strike K. The 
associated conditional density process {ft(x)} for the random variable At is defined by 

[" f t (x)te = EP[l{y> A T }\JF t ). (1.3) 

J — oo 

Then for the asset price we have 

A t = / xf t (x)dx, (1.4) 



and the corresponding option prices are given by 

C tT (K) = f(x- K) + f t (x)dx. (1.5) 



Instead of modelling {A t } and then determining {C t T{K)}, our strategy is to model the 
conditional density process. Then the underlying asset price process and the associated 
option prices are determined by (1.4) and (1.5). Roughly speaking, the idea is to model 
{ft{x)} in such a way that it contains some parametric freedom that can be calibrated to a 
specified range of initial option prices. Models for conditional densities have been considered 
in various contexts in finance. These include for example applications to interest rates (Brody 
& Hughston 2001a,b, 2002, Filipovic et al. 2010), and to credit risk (El Karoui et al. 2010). 
Although mostly different from what has previously appeared in the literature, our approach 
to modelling the volatility surface is similar in spirit in some respects to that of Davis (2004). 

Let us consider in more detail the class of assets that will form the basis of our inves- 
tigation. We introduce a probability space (f2, J 7 , P) with filtration {J-"t}t>o, where IP is the 
"physical" measure and {Ft} is the market filtration. We assume that price processes are 
adapted to {Ft}- We assume the absence of arbitrage, and the existence of an established 
pricing kernel {lit} associated with some choice of base currency as numeraire. We work in 
the setting of a multi-asset market, and do not assume that the market is complete. 

We write {^}i=o,i,...,Ar for the price processes of a collection of non-dividend-paying trad- 
able financial assets. Prices are expressed in units of the base currency. We refer to asset i 
as A\ We model the {A\} as Ito processes, and for each A 1 we require that 

n s A l s =E ¥ [ir t Al\F s ] (1.6) 

for s < t. Such an asset is characterized by its value A T at some terminal date T. In some 
situations it is useful to regard the asset as offering a single payment at T. In that case 
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{Al\ < t<T represents the price process of the asset that offers such a payment. In other 
situations we can consider At as being a "snapshot" of the value of the asset at time T. 
Usually the context will make it clear which meaning is intended. 

Let {A® } be a money-market account in the base currency, initialised to unity. If we 
set p t = n t At for t > 0, it follows that {p t } is a P-martingale. A standard argument shows 
that {p t } can be used to make a change of measure. The resulting measure Q° is the "risk- 
neutral" measure associated with the base currency, and has the property that if the price 
of any non-dividend-paying asset is expressed in units of the money-market account, then 
the resulting process is a Q°-martingale. Thus, for each i we have 



A] 



A° s E ( 



A° t 



;i.t) 



A similar situation arises with other choices of numeraire. Specifically, for any non- 
dividend-paying asset A % of limited liability (A\ > 0), with price {A l t }, there is an associated 
measure Q* with the property that if the price of any non-dividend-paying asset is expressed 
in units of A 1 then the result is a Q J -martingale. Thus for < s < t and for all i, j for which 
the price of A 1 is strictly positive we have: 



At = AlK- 



'4 






At 



(1. 



Bearing these points in mind, we observe that the option pricing problem can be formu- 
lated in the following context. We consider European-style options of a "Margrabe" type, 
for a pair of non-dividend-paying assets A % and A\ where the option-holder has the right at 
time t to exchange K units of A 1 for one unit of AK The payoff H\ 3 of such an option, in 
units of the base currency, is of the form 



HnK) = {A\-KA\) + . 

The value of the option at s < t, expressed in units of the base currency, is given by 

1 



Cll(K) 



7T, 



ix t {Ai-KA\y \F S 



;i.9) 



;i.io) 



If A % is of limited liability, then the option value, expressed in units of A 1 , is a Q*-martingale: 

eg (JO 



Ai 



Ai 



-K 



1.11 



This relation can be expressed more compactly as follows. Write A s for the price at time 
s < t of a "generic asset" expressed in units of a "generic numeraire", and C st {K) for the 
price, in units of the chosen numeraire, at time s < t, of a t-maturity .fT-strike option. Then 
the option payoff is given (in numeraire units) by 



H t {K) = {A t - K)~ 



(1.12) 
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The value of the option at time s < t is 

C st (K)=¥®[(A t -K) + \J?.], (1.13) 

where Q is the martingale measure associated with the numeraire. By "generic" we mean 
any choice of a non-dividend-paying assets A and A 1 such that A is of limited liability. 

Standard options are not included in the category discussed above. A standard call has 
the payoff H t = (A t — K) + where A t is the price of the underlying at time t in currency 
units, and K is a fixed strike in currency units. This is not an option to exchange K units of 
a non-dividend-paying asset for one unit of another non-dividend-paying asset. One should 
think of the fixed strike as being K units of a unit floating rate note. The only asset that 
maintains a constant value (in base-currency units) is a floating rate note — and such an asset 
pays a dividend. The dividend is the interest rate. In a given currency, by a floating rate 
note we mean an idealised note that pays interest continuously (rather than in lumps). The 
associated dividend is the short rate. One might argue that the strike on a standard option 
is cash, and that cash is a non-dividend-paying asset: this point of view leads to paradoxes. 
In the standard theory we regard cash as paying an implicit dividend, a convenience yield, 
in the form of a liquidity benefit equivalent to the interest rate. In summary, a standard 
option is an option to exchange the asset with a floating-rate note, which pays a "dividend" . 

Thus a standard option is a complicated entity — it is an option to exchange a certain 
number of units of a dividend-paying asset for one unit of a non-dividend-paying asset. It 
is more logical first to examine an option based on a pair of non-dividend-paying assets. In 
the literature this approach is implicitly adopted through the device of "setting the interest 
rate equal to zero". In that situation the floating rate note is non-dividend-paying; thus, 
the setting we operate within includes the zero-interest case. One would like to tackle the 
general problem of an option to exchange K units of one dividend-paying asset for one 
unit of another dividend-paying asset (a "standard" foreign exchange option falls into that 
category); but, unless the dividend (or the interest rate) systems are deterministic, this is a 
more difficult problem than the one we propose to consider here. 

The structure of the paper is as follows. In Section 2 we derive a dynamical equation 
for the conditional density, which we call the "master equation", given by ( I2.13p . An in- 
tegral form of the equation is presented in (I2.14p . which forms the basis of the solutions 
presented in later sections. In Section 3 we specify the general form we require the volatility 
structure of the conditional density to take, and give a characterization of what we mean 
by a "conditional density model" for asset pricing. In Section 4 we consider in detail the 
class of models for which the volatility structure of the conditional density is a deterministic 
function of two variables. This family of models admits a complete solution by use of a 
filtering technique. The resulting asset prices exhibit a stochastic volatility that is adapted 
to the market filtration but that is not in general of the local- volatility type. In Section 5 we 
consider the case when the volatility structure is linear in the terminal value of the asset. In 
that case the resulting models are Markovian, and can be calibrated to an arbitrary initial 
density. In Section 6 an alternative representation of the semi-linear case is presented using 
a Brownian-bridge technique. In Section 7 we show that the Bachelier model and the geo- 
metric Brownian motion model arise as special cases of the semilinear models, for particular 
choices of the initial density. We conclude in Section 8 with the calculation of option prices. 
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2 Conditional density processes 



The market is understood as having the setup described in the previous section. We have 
a probability space (Q, J 7 , Q) with filtration {J^}. A non-dividend-paying limited-liability 
asset is chosen as numeraire, and all prices are expressed in units of that numeraire. The 
measure Q has the property that the price process of a non-dividend-paying asset, when 
expressed in units of the numeraire, is a martingale. We refer to Q as the martingale 
measure associated with this numeraire. 

We fix T > 0, and assume for < t < T the existence of an ^-conditional Q-density 
ft{x) for At- In our applications we have in mind the cases iel and x £ R + , but it is useful 
to be flexible as regards the choice of the domain of the density function. In what follows we 
treat the case x £ M, and leave it to the reader to supply the necessary adjustments for other 
domains. Thus we assume the existence of a density process {f t (x)}, x £ R, < t < T, such 
that (11.31) holds. It follows that for any bounded, measurable function g(x), x £ K, we have 



/ g(x)f t (x)dx = E[g(A T )\^]. (2.1) 
Jm 



We ask that {ft(x)} should have those properties that follow heuristically as consequence of 
the formula obtained by formally differentiating (II. 3p . namely: 

f t (x)=E[5(x-A T )\F t ]. (2.2) 

In particular, we require the following: (a) that for each x £ K. the process {ft{x)}o<t<T is 
an {J^j-martingale, and hence f s (x) = E [ft{x) | J-" s ] for < s < t < T; (b) that the price of 
the asset can be expressed in terms of the density by (II. 4p for t < T; and (c) that 

lim / xf t (x)dx = A T . (2.3) 

We also assume, where required, that expressions analogous to (11.41) can be written for 
claims based on A T . For example, if C t T(K) denotes the price at t of a T-maturity, i^-strike 
European call option, then we assume that (II. 5p holds. 

In the applications that follow, we introduce a Q-Brownian motion {Wt}, which we take 
to be adapted to {J-'t}, and specialize to the case for which the dynamical equation of {ft(x)} 
is of the form t 

ft{x) = f (x) + / <j{{x) f s (x)dW s , (2.4) 



for some process {a{(x)}, x £ 1R, < t < T, representing the volatility of the density. It 
follows that t 

A t = A + [ a^dW s , (2.5) 



n 



where A = J R xf (x)dx, and where the volatility of A t is given by 

of = [ xa{(x)f t (x)dx. (2.6) 



For simplicity, we consider in this paper the case where {Wt} is a one-dimensional Brownian 
motion. The extension to the multi-dimensional situation is unproblematic. 
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Lemma 2.1. The normalisation condition 



ft(x)dx = 1 



holds for all t G [0,T) if and only if there exists a process {a t (x)} such that 



(2.7) 



a{(x)=at(x) / f t (y)dy- / a t (y)f t (y)dy 



(2.1 



for almost all t G [0, T), and we have the initial condition 



/o(x)cb = 1. 



Proof. First we show that (12. 8p and (12.91) imply (12. 7p . Starting with 
with respect to x to obtain 



(2.9) 



we integrate 



f t (x)dx= / f (x)dx+ / / ni(,-)./;(,-)dr(lli;. 



(2.10) 



Inserting (I2.8P and using (I2.9p we obtain (12. 7p . Conversely, if we assume (12. 7p then (12. 9p 
holds as well, and hence 



a{{x)f s {x)dxdW s = 0, 



(2.11) 



from which it follows that 



a{{x)f t {x)dx = (2.12) 
for almost all t G [0,T), and thus that (12.81) holds for almost all t G [0, T) for some <7t(x). □ 

Thus, once we specify {a t (x)} and {/o(^)}, the dynamical equation for the density — the 
so-called master equation — takes the form 



ft(x) = f (x) + 



t r 



a s (x) - / (T s (y)f s (y)dy 



fs{x) dW s 



(2.13) 



Lemma 2.2. The conditional desity process {ft(x)} satisfies the master equation (I2.13P 
with initial density {fo(x)} and volatility structure {&t(x)} if and only if 



fo(x) exp ( f* a s (x)dZ s - \ f* a s 2 (x)ds 
/ R fo(y) exp (/„* a s (y)dZ s - \ f* a s 2 (y)ds) dy 



where 



- / / <r.(y)fs(y)dyds. 



(2.14) 



(2.15) 



(3 



Proof. Writing ( I2.13P in differential form, we have 

df t (x) = f t (x) (a t (x)-(a t ))dW t 
where for convenience we write 

(<?t) = / ar t {x)f t (x)dx 



for the conditional "mean" of the volatility. We integrate (I2.16P to obtain 



ft(x) = f (x)exp 



(a s (x) - (a s )) dW s -\ (a s (x) - (a s )) 2 ds 



Expanding the exponent, we have 



ft(x) = f (x) 



cxp 



£ a s (x) (dW s + (a s )ds) - \ jj a s \x)ds 



exp 



f (vs) (dW. + (a.)ds)-\fi(a.)*ds 



Then we introduce a process {Z t } by writing 

Z t = W t + [ (a s )ds, 
Jo 



and it follows that 



ft{x) = f (x)- 



exp 


Io^s(x)dZ s - \f Q a s 2 (x)d 


s 


exp 







Applying the normalization condition (12. 7p to equation f)2.2ip above, we see that 



(a s }dZ s - | / (a s ) 2 ds 



(2.16) 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



(2.21) 



/ f (x)exp( [ o s (x)dZ s -\ [ a s 2 {x)ds) dx. (2.22) 
Jr \Jo Jo / 



It follows that (I2.2ip reduces to (I2.14p . Conversely, if ft{x) is given by (I2.14p . then it is a 
straightforward exercise in Ito calculus to check that the master equation is satisfied. □ 

3 Conditional density models 

We are in a position now to say more precisely what we mean by a "conditional density 
model" . In doing so, we are motivated in part by advances in the study of infinite-dimensional 
stochastic differential equations. By a "model" for the density process we understand the 
following. We consider solutions of the master equation ( I2.13P satisfying the normalization 
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condition (12 .7p . in conjunction with the specification of: (a) the initial density fo{x); and 
(b) the volatility structure {(J t (x)} in the form of a functional 



a t (x) = ^[f t (-),t,x]. (3.1) 

For each x and t the volatility at(x) depends on ft(y) for all y £ R. Hence (12.131) . thus 
specified, determines the dynamics of an infinite-dimensional Markov process. 

The initial density fo(x) can be determined if one supplies initial option price data for 
the maturity date T and for all strikes K e R. In particular, we have 

Cot(K) = E[(A T - K) + ] = f(x - K) + f (x)dx. (3.2) 

Jr 

By use of the idea of Breeden & Litzenberger (1978) we see, in the present context, that for 
each value of x one has 

/oW = — ^ — • (3-3) 

Here fo(x) is not generally the risk-neutral density, but rather the Q-density of the value at 
time T of the asset in units of the chosen numeraire. For simple practical applications, we 
take — in common with much of the literature — the numeraire to the be the money-market 
account, and set the interest rate to zero. Then fo(x) is the risk-neutral density. Once fo(x) 
has been supplied, the choice of the functional $ determines the model for the conditional 
density: we give some examples later in the paper. 

In practice, one would like to specify <3> modulo sufficient parametric freedom to allow the 
input of additional option price data. What form this additional data might take depends on 
the nature of the market and the class of valuation problems being pursued. For example, 
a standard problem would be to look at a limited-liability asset and consider additional 
data in the form of initial option prices for all strikes in IR + and all maturities in the strip 
< t < T. We require that $ should be specified in such a way that once the data are 
provided, then $ is determined and the "master equation" provides an evolution of the 
conditional density. Once we have the conditional density process, we can work out the 
evolution of the option price system for the specified strip, and hence the evolution of the 
associated implied volatility surface. 

The data do not have to be presented exactly in the way specified in the previous 
paragraph — there may be situations where more data are available (e.g., in the form of 
barrier option prices or other derivative prices) or where less data are available (less well- 
developed markets). One should think of the parametric form of $ as being adapted in a 
flexible way to the nature of a specific problem. The philosophy is that there are many 
different markets for options, and one needs a methodology that can accommodate these 
with reasonable generality. 

It should be evident that an arbitrary solution to the master equation need not be a 
density process for an JVmeasurable random variable — additional assumptions are required 
concerning the nature of the volatility structure in order to ensure that ft{x) converges in an 
appropriate sense to a suitable Dirac distribution. In the examples given, we indicate how 
this can be achieved, in various situations, by the choice of the volatility structure. 
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4 Models with deterministic volatility structures 



We proceed to present a rather general class of conditional density models characterized by 
a deterministic volatility structure. These are models for which {a t (x)} is of the form 



(TAX) 



v(t, x) 



(4.1) 



for some deterministic function v(t, x) defined for appropriate values of t and x. We are able 
to give a more or less complete construction of such models in the form of a "weak" solution 
of the master equation. By a weak solution we mean that on a probability space (Q, J 7 , Q) 
we construct a filtration {J"t}, a Brownian motion {W^}, and a conditional density process 
{ft{x)} satisfying the master equation and having the desired properties. More specifically, 
let T e (0, oo) be fixed, and let fo(x) : K — > M + be a prescribed initial density. Let v(t, x) 
be a function on [0,T) x M satisfying 



(4.2) 



v (s, x) 2 ds < oo 

for t < T and x6R, and let ^(t) be a function on [0, T) such that 

lim7(t) = 0, lim7(t) / v (s, x)ds = g(x), 

J 

for some invertible function g(x) on IL We have the following: 



(4.3) 



Proposition 4.1. Let X have density fo(x), and let {B t } be an independent Brownian 
motion. Let {J 7 *} be the filtration generated by the "information" process {It} defined by 



h = B t 



v(s, X)ds. 



(4.4) 



Let {ft(x)} be defined by 
ft(x) 

and define {W t } by setting 



fo(x)exp 


Jl ) v(s,x)dl s - 


lf*v 2 (s,x)ds 




SZo My) exp 


J*v(s,y)dl s - 




dy 



(4.5) 



W t = I t 



[ E Q [v(s,X) | JF S ] ds. 
Jo 



(4.6) 



Then: (a) the random variable X is J^-measurable, and ft(x) is the associated conditional 
density; (b) the process {W t } is an {J^j-adapted Brownian motion; and (c) for t G [0, T), 
the density process {ft{x)} satisfies the master equation 



ft(x) = f (x) + / f a (x) 



v(s,x) - / v(s,y)f s (y)dy 



(4.7) 
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Proof. It follows from (14. 3 p that g(X) = \im t ^ T ^{t)I t is J^-measurable. Since g(x) is 
invertible, we conclude that X is J^-measurable. Now let the filtration {Gt} be defined by 



Qt — cr ({B s } < s < t ,X) . 

Clearly Tt C Gt- The random variable X is (^-measurable, and {B t } is a ({Gt}, 
motion. We introduce a ({Gt}, Q) -martingale {M t } by setting 



M t = exp 



f v(s,X)dB s -± [ v 2 (s,X)dt 
Jo Jo 



(4.8) 
-Brownian 

(4.9) 



for t G [0,T), and we let the probability measure B be defined by 

dB 

M t . 



(4.10) 



We observe that dl t = dB t + v(t,X)dt. By Girsanov's theorem, {I t } is a {{Gt}, B)-Brownian 
motion. We note that {M t -1 } is a ({&}, B)-martingale. Let if be a bounded measurable 
function on R. Since H{X) is ^-measurable, we have the generalised Bayes formula 



[H(X)\F t \ 



E B [Mt X H{X) | Ft] 
E B [M t _1 1 F t ] 



(4.11) 



Let us work out the right-hand side of this equation. To this end we show that X and /( are 
B-independent for all t. In particular, we show that the generating function 



factorises. We have: 
E B [exp (yl t - 



E B [exp (yl t + zX)] 



zX)]=E Q [M t exp(yI t + zX)} 



E { 



exp 



v(s,X)dB s 



v 2 (s,X)ds 



+y 



B, 



v(s, X)ds 



exp (zX) 



E Q [mt exp (\y 2 t) exp (zX)] 



where 



m t = exp 



-v(s,X) + y} 2 ds 



-v(s,X)+y]dB. 

By use of the tower property and the independence of {B t } and X under 
E B [exp (yl t + zX)] 



(4.12) 
(4.13) 

(4.14) 
(4.15) 

(4.16) 



we have 



E Q [m t exp (\y 2 t) exp (zX)] , (4.17) 
exp (\y 2 t) E Q [E Q [m t exp (zX) \ X}] , (4.18) 
exp (\y 2 t) E Q [E Q [m t \ X] E Q [exp (zX)]] . (4.19) 
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One observes that E^[m t | X] = 1. Thus we obtain the desired factorization: 

E B [exp (yl t + zX)} = exp (\yH) E Q [exp (zX)} . (4.20) 

Now that we have shown that X is IB-independent of I t (and thus of Ft), we can work out 
the right-hand side of ( 14. lip . We have: 



Or , x, , E \Mf H(X) Ft\ 
E Q # (X) LFJ = ; y ' tJ 

e» Af t - j-j 



E B 


H(X) exp ( J * v(s, X)dl s - \ W 2 ( S , X)ds 






E B 


exp (/ ^( S ,X)dJ s - l/ V(s,X)ds) 


I 







/!L fo{x)H(x) exp ( J r v(s, x)dl s - \ f* v 2 (s, x)ds ) dx 



IZo My) ex P ( Jo ^) dJ ^ - \ Jo ^ 2 ( s ' ?/) ds ) d y 



(4.21) 
(4.22) 

(4.23) 



In particular, setting H(X) = 1(X < x), we deduce that ft(x) is the J-" t -conditional density 
of X, as required. That proves the first part of the proposition. 

Next we show that {Wt} is an ({J-t}, Q)-Brownian motion. We need to show (1) that 
(dW t ) 2 = dt, and (2) that E Q [W U | F t ] = W t for < t < u. The first condition is evidently 
satisfied. The second condition can be shown to be satisfied as follows. For simplicity, we 
suppress the superscript Q. We have 



E[W^ u |7i]=E[/ tt |J=i]-E / E[v(a,X)\ F s ]ds\ F t 

Jo 

First we work out E [I u \ Ft]. Since {B t } is a ({Gt}, Q)-Brownian motion, we have 



(4.24) 



E[I U \ Ft] 



E 



B u + v(s,X)ds\F t 
Jo 



E [B u | F t ] + E 
E [E [B u | &] | Ft] + E 
E [S t | F t ] + E 



/ u(s,X)ds| 7* f 

/ ?;( S ,X)d S |.F t 



f (s, X)ds I 



(4.25) 
(4.26) 
(4.27) 
(4.28) 



We insert this intermediate result in (I4.24p to obtain 



E[W U \F t ]=E [I u \F t ]-E 
= E[B t \ Ft] + E 



/ E[v(s,X)\F s ]ds\F t 
Jo 



\F t 


-E 











/ E[v{s,X)\F s ]ds\F t 
Jo 



(4.29) 
(4.30) 
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Next we split the integrals in the last two expectations by writing, 



E[W U | Ji] = E [B t | F t ] + E 



/ v(s,X)ds\F t 
Jo 



E 



/ EK S ,JQ| F s ]ds\ F t 
Jo 



E 



v(s,X)ds\F t 
E[v(s,X)\F s ] ds\F t 



Observing that 



and that 



E 



E [B t | F t ] + E 



v(s, X)ds I J-i 



E 



v(s, X)ds I J-i 



E[v(s,X)\F s ] ds\F t 



(4.31) 

(4.32) 
(4.33) 



we see that the expectation E[WJ Ft] reduces to 



E[W U | Ft] = h - E 



/ E [v(s, X) | F s ] ds\F t ]=It- [ E X) | F s ] ds = W t . (4.34) 
Jo J </o 



That shows that {W t } is an {J-" 4 }-Brownian motion. An application of Ito calculus shows 
that the density process (14. 5 j) satisfies the SDE (14.71) . □ 



Remark 4.1. For simplicity we have presented Proposition 14.11 for the case of a one- 
dimensional state space. It should be clear from the proof how the results carry over to 
conditional densities on higher-dimensional state spaces. 

Remark 4.2. For simulations of the dynamics of the conditional density process, the fol- 
lowing alternative representation for (14 . 5 1) may prove useful: 



fo(x) exp (/ [v(s, x) - v(s, X)] dB s -\J Q [v(s, x) - v(s, X)] ds 
/ R /o(y)exp (j*[v(s,y) - v(s,X)]dB a - |/* [v(s,y) - v(s, X)} 2 ds) dy 



(4.35) 



The unnormalised density — the numerator in (I4.35P — is for all x £ R conditionally log- 
normal given X. The simulation of the density requires only the numerical implementation 
of the standard Brownian motion and of the random variable X. 



Remark 4.3. The density models with deterministic volatility structure presented in Propo- 
sition 14.11 can be extended to a class of models that satisfy the following system: 



Let fo(x) : R — > R + be a density function. A filtered probability space (f2, F, {F t },Q) 
can be constructed along with (i) an J^-measurable random variable X with density fo(x), 
(ii) an {J-j}-adapted density process {ft(x)}, and (iii) an {J-i}-adapted Brownian motion 
{Wt}, such that (a) for some function j(t) on [0, oo), (b) for some function g(x) that is 
invertible onto R, and (c) for some suitably integrable function v(t,x) on [0, oo) x R with 
the properties 

lim y/t^it) — 0, lim 7(i) / v(s,x)ds = g(x), (4.36) 
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the following relations hold for all t G [0, oo). We have Q [X e dx \ J- t ] = ft{x) dx and 



f t (x) = f (x) + / f a (x) 







v(s,x) - / v(s,y)f s (y)dy 



dW s . (4.37) 



Remark 4.4. The construction of conditional density models admits an interpretation as 
a kind of a filtering problem. The process (I4.4p plays the role of an "observation process" , 
and (14.61) has the interpretation of being an "innovation process". See [T], fTS], |16j . 

Remark 4.5. It is reasonable on a heuristic basis to expect that the general deterministic 
volatility structure model can be calibrated to the specification of an essentially arbitrary 
volatility surface. In particular, the parametric freedom implicit in a deterministic volatility 
structure coincides with that of a volatility surface. The situation is rather similar to that 
of the relation arising in the Dupire (1994) model between the local volatility (which is 
determined by a deterministic function of two variables, one with the dimensionality of 
price and the other with that of time) and the initial volatility surface (which represents 
a two-parameter family of option prices, labelled by strike and maturity). The precise 
characterisation of such relations constitutes a non-trivial and important inverse problem. 

Remark 4.6. In the general deterministic volatility structure model, the dynamics of the 
underlying asset price {^4*} are of the form dA t = V t dW t , with 

V t = xv(t,x)f t (x)dx - / xf t (x)dx / v(t, x)f t (x)dx, (4.38) 
Jr Jr Jr 

where ft{x) is given by (14.51) . Thus the absolute volatility Vt at time t takes the form of a 
conditional covariance between X and v(t, X). In general, {V t } is an {J^j-adapted stochastic 
volatility process that cannot be expressed in the form T,(t,A t ) for some function £(£, x), 
and the dynamics do not constitute a simple diffusion of the Dupire (local volatility) type. 
In the "semilinear" case v(t,x) = aTx/(T — t), however, the associated information process 
is Markovian, and V t can indeed be expressed in the form E(t, A t ). 



5 Semilinear volatility structure 

We consider in this section the case where At has a prescribed unconditional density /ot(^), 
and construct a family of conditional density processes {ftr} that solve the master equation 
(I2.13P over the time interval [0, T). The filtration with respect to which {ftr} is defined will 
be constructed as follows. We introduce a process {£,tT}o<t<T given by ^ t r = cr A T t + (3 tT , 
where a is a constant and {/3tr}o<t<r is a standard Brownian bridge, taken to be independent 
from At- We assume that {Ft} is given by 

^^(U<k| (5-1) 

Clearly {C,tr} is {J"t}-adapted, and At is J-r-measurable. It is shown in Brody et al. (2007, 
2008) that {^t} is an {J-" t }-Markov process (see also Rutkowski & Yu 2007). 
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Proposition 5.1. Let the initial density Jot(x) be prescribed, and let the volatility structure 
be of the semilinear form 

T 



a tT {x) = a 



T-t 



x. 



(5.2) 



for < t < T, and let {T t } be defined by (15. ip . Then the process {W t }o<t<T defined by 



{(TTE[A T \U]-U) ds 



! T-s 

is an {J^j-Brownian motion, and the process {/^(x)}, given by 

/otO) exp [f^ (a£ tT x - \ a 2 x 2 t)] 



(5.3) 



ftT(x) 



IZo My) exp [f^ (a£ tT y - \ a 2 y 2 t)] dy'' 
satisfies the master equation (I2.13P with the given initial condition. 



(5.4) 



Proof. The fact that {W t }o<t<T is an {J^j-Brownian motion is shown in Brody et al. (2007, 
2008). We work out K[A T \ £ tT ] by use of the Bayes formula, 



where 



ftr(x 



IE [At I Ctr] = [ x f tT (x)dx, 
Jr 

hT{x)p(£ t T I A T = x 



(5.5) 



(5.6) 



f M for{y) P (&r I A T = y) dy' 

Here p(£tT | At = x) is the conditional density of £tT given the value of At- We observe that 
£tr is conditionally Gaussian: 



P (&r \A T = x) 



T 



2irt(T -t) 
Thus the density process is given by 

ftAx) 



exp 



T 



2 t(T - 1) 



(5.7) 



f 0T (x)exp 


2 t(T-t) atX ) \ 




IZo My) ex P 


I 2 t(T- t ) (fa at y) 


dy 



(5.1 



The last expression can be simplified after some rearrangement so as to take the form 

_ f 0T (x) exp [ggj (ajtTX - \ a 2 x 2 t)] 

* JZc My) ex P [t^i ( a &Ty - \ ° 2 y 2t )] &y ' 

With this result at hand, we can write the process {W t }o<t<T i n the form 



(5.9) 



W t = i tT - 



T-s 



crT x f sT (x)dx - £ sT ds. 



(5.10) 
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We recall that the master equation (I2.13P can be written as (I2.14p . We shall prove that (15 .4p 
satisfies (12 . 1 3[) by showing that (I2.14p reduces to (15. 4p if we insert (15. 3p in (I2.20p and choose 
the volatility structure to be (15. 2p . For the process {Z t } in (I2.20p we obtain 



T^~s dS ' 



(5.H) 



The next step is to insert {Z t } in the exponent 

/ a sT (x)dZ s - \ [ a sT (x) 2 ds (5.12) 
Jo Jo 

appearing in equation (I2.14p . Expression (I5.12p can be simplified by use of (15. 2p to give 



^r- — ( d£ sT 1 



T-s 



T-s 



To derive this result, we make use of the relation 



sT 



, T-s T-t J (T-s)" 
With equation f[5TT3"j) at hand, we see that fl2TT4l) reduces to f[5T4l if fl5T2l holds. 



T 

= 7jT—^ {vxitT - \ cr 2 x 2 t) . 

(5.13) 

(5.14) 
□ 



ds. 



6 Semilinear volatility: Brownian motion approach 

We proceed to show how the models constructed in Section are related to the density 
models with deterministic volatility structure treated in Section HI In particular, we consider 
a deterministic semilinear volatility function of the form 

T 

v(t, x) = a 



T-t 



x, 



where < t < T. For this volatility function the process {It} has the dynamics 



dl t = a 



T 



T-t 



Xdt + dB t . 



We are thus able to work out the exponent 

/ v(s,x)dl s — | / v 2 (s,x)ds 
Jo Jo 

in equation (14 .5p making use of (16. ip and (16. 2p . We have: 



v(s,x)dl s — \ J v 2 (s,x)ds 



T 



T-t 



ax 



(T-t) 



f^ + T(T-t)aX f 
Jo T ~ s Jo 



(6.1) 



(6.2) 



(6.3) 



d.s 



(T-s) 



1t>2 2 2 

■kl a x 



ds 



(T-sf 

(6.4) 



15 



The first integral gives rise to a ({Gt\, Q)-Brownian bridge {At} over the interval [0,T]. 
More specifically, we have 

(3tT = (T-t) f^. (6.5) 



(6.6) 



o T 

The deterministic integral in (16. 4p gives 

"* ds t 
(T - s) 2 ~ T(T — t) ' 

Armed with these results, one can write ( 16. 3 p as follows: 

J v(s,x)dl a -±j v\s,x)ds= ^—ax(aXt + f3 tT )-^j?— t (T 2 x 2 t. (6.7) 

Let {£tr} be defined for t G [0,T] by £ t T = cr Xt + f3 t r- Then for (16. 7p we obtain 

/ v(s, x)dl s - \ I v 2 (s, x)ds = T (crx^tT - \ cr 2 x 2 t) . (6.8) 
Jo Jo 1 — t 

We conclude that the conditional density process {ft(x)} in (14 .5p reduces to the following 
expression in the case for which the volatility structure is given by (16. ip : 

fo (x) exp \t3—; (cr x £ t r — h cr 2 x 2 t] 1 
f t (x) = ^ v 1 lT - f v 2 11 . (6.9) 

IZo fo(v) ex P { a v£tT - \ cr 2 y 2 t)] dy 

From equation (16. 8 p we see that {^t} takes the role of the information process that generates 
{J-f}. Since {/3*t} vanishes for t = T, the random variable X is "revealed" at T. Thus X is 
J-^-measurable, and {£tr} is the process generating the information-based models of Brody 
et al. (2007, 2008). This conclusion is supported by the following construction. 

We consider the measure B defined in (14.101) . Under B the process {I t } is a Brownian 
motion over the interval [0, T). We construct a B-Brownian bridge by use of the B-Brownian 
motion {I t } as follows. On [0, T) we set 

itT = (T-t) /" -J_dJ s . (6.10) 
Jo 1 ~ s 

Next we recall definition (I4.4p and insert this in the expression above. The result is 

itT = (T-t) f ^ + {T-t) f -J—v(s,X)ds. (6.11) 
Jo T-s Jo T-s 

The first integral defines a {{Gt}, Q)-Brownian bridge over the interval [0, T) which we denote 
{fitr}- For the volatility function we set 

v(t,x) = a 7 ^— t x. (6.12) 

This leads to ^ 

£ tT = aX(T- t)T [ 1 ds + At, (6.13) 
Jo ( T - s ) 

and thence to (6.8). 
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7 Bachelier model 



The Bachelier model is obtained by setting A t = "fWt where 7 is a constant. We shall show 
that the class of models defined by Proposition 15 . II contains the Bachelier model. We consider 
a random variable At associated with a fixed date T. We assume that At ~ N[Q, l/(Ta 2 )}, 
where N[m,v] is the class of Gaussian random variables with mean m and variance v. In 
the notation of Section [5j we have 



a 



exp ( — ~ cr 2 Tx 2 ) . 



(7.1) 



We recall the process {^t}o<j<t defined by £ t T = cr A T t + f3 tT . If A T ~ -/V(0, l/(Tcr 2 )), then 
{£tr} is an {J-^j-Brownian motion over [0, T]. This is because {£tr} is a continuous Gaussian 
process with £ t = and Cov[£ s t, £tr] = s for < s < t < T. We recall the definition of the 
Brownian motion {W t } associated with {^t}, given by (15. 3p . Since for At ~ ^(0, 1/Tcx 2 ) 
the process {£tr} is a Brownian motion, it follows that 



E [A T I £ sT ] 



Thus we see that Wt = £tT- As a consequence we have 



A t = E [A T I &t] 



±E[Wr|Wi] = iwi. 



(7.2) 



(7.3) 



Hence to match the Bachelier model with an element in the class of models constructed in 
Section [51 it suffices to set cr = 1/(7 T). 

Proposition 7.1. The conditional density process {f^{x)} of the Bachelier price process, 
defined over the interval [0, T), given by 



exp 






TOO 

J-oo eX P 


-It^M^-tW,) 2 


dy 



(7.4) 



is a special case of the family of the models of Proposition 15. 1\ and is obtained by setting 



/or (a?) 



exp (— ~ cr 2 Tx 2 ) , cr tT (x) = a 



T 



;r. 



and 



Proof. We insert (17. 5p in (15. 4p . Completion of squares gives 

/or(ic) exp (cr^ T a; - | erVt)] 



V(7T). 



(7.5) 



/^o /or(y) exp (cr^Ti/ - I a 2 ?/ 2 t)] cfy 



exp 


l^T 2 ( 1 t \2] 




poo 

J-oo eX P 


1ct2T2 f« 1 ^^ 2 


dy 



(7.6) 



Recalling that ^ t T = W t , and setting a = 1/(7 T), we obtain the desired result. □ 
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Remark 7.1. Let the initial density for{x), the volatility structure cr tT (x), and the param- 
eter cr be given as in (17.51) . Then the Bachelier conditional density {f^r(x)} satisfies the 
master equation (12.131) . and {W t } coincides with {^t}- 

Remark 7.2. Suppose we chose an asset price model with a certain law. Then we know 
that we can derive the corresponding conditional density process where the related volatility 
structure and initial density are specified. We may then wonder how the conditional density 
transforms, and what the new volatility structure looks like, if we consider a new law for 
the asset price model. For instance, we may begin with the Bachelier model and ask what 
is the conditional density and volatility structure associated with a log-normal model. We 
present a "transformation formula" for the conditional density. This result allows for the 
construction of a variety of conditional density processes from a given one. Let {ft{x)} solve 
(I2.13p . and let ip : R — > M be a C 1 -bijection. Then it is known that if X has conditional 
density ft(x) then Z = ip(X) has conditional density gt(z) given by 

9 ' (2) = WSFWY (77) 

where ip~ 1 (z) is the inverse function and ip'(x) is the derivative of ip(x). The conditional 
density {g t (z)} satisfies 

dg t {z)=g t {z)[v{t,z)-{v t )]dW t , (7.8) 
where u(t, z) = v{t,ip~ l {z)) and 

(u t )= [ v(t,y)g t (y)dy. (7.9) 

We see that the volatility structure v(t,x) associated with ft(x) becomes the volatility 
structure v(t, z) associated with gt(z). For example, consider the Bachelier model where 
A T ~ N(0, 1/(T<7 2 )). The associated initial density and volatility structure are given in 
(17. 5p . Now suppose that Z = exp(A T ), so ip(x) = exp(x). The price process {A t } is then 
given by the log-normal model 

A t = exp( T W t ), (7.10) 

where W t = £tT- It follows by (17. 7\i that the conditional density process {gtr(x)} associated 
with the log-normal price process (I7.10P is 



gtT{z) 



exp 






r oo 

z Jo ex P 


-l^fc) (My) - Wtf 


dy 



(7.11) 



for z > 0. Indeed we see that gtr{z) is the log- normal conditional density. The associated 
volatility structure is 

u(t, z) = v(t, \n(z)) = a tT (z) = a ^-^ \n{z). (7.12) 
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8 Option prices 



We consider a European-style call option with maturity t, strike K, and price as determined 
by equation (j!.13p . The price process {A t } < t <oo of the underlying asset is given by 



-4, 



E Q [X | F t ] 



JZo X fo( X ) eX P (/o X ) dI * ~ \ Jo t;2 ( S ' X ) ds ) dx 



IZc fo(y) ex p ( So y) dI * - \ Jo ^ 2 ( s ' y) ds ) d y 



i.r 



where {^t}o<s<t<oo is generated by (14. 4p . We recall that {B t } is a 
motion, and that there exists an ({-7^}, Q)-Brownian motion {W t } such that 



-Brownian 



dl t = dW t 



(v t )dt, 



(8.2) 



where the bracket notation is defined by (12.171) . We introduce a positive {{J~"t}, Q)-martingale 



A t = exp / (v s )dW s + \ \ (v s ) 2 ds 
\Jo JO / 

which induces a change of measure from Q to measure Q* given by 

= At- 



Ft 

The Q*-measure is characterised by the fact that {I t } is an ({J^}, 
We observe that by the relationship 

/o(^)exp( / v(s,x)dl s — \ I v 2 (s,x)ds\dx 



(8.3) 
(8.4) 

-Brownian motion. 



Jo 
exp 

:A,, 



(v s )dl s 



s 2 

o Jo 



(v s ) 2 ds 



(8.5) 
(8.6) 



we can use the denominator in (18. ip to write the option price C s t in terms of a conditional 
expectation taken with respect to Q* under which {I t } is a Brownian motion. Equation 
(18.61) is obtained by applying the relationship (18. 2p . We then have 



{N t A^-K) + 
A" 1 (N t - KA t 



E 1 



J 7 , 



where 



N t 
A t 



AJ 1 VP [(N t -KA t ) + |J-J, 



x fo(x) exp v(s,x)dl s — \ J v 2 (s,x)ds*j dx, 
fo(x)exp I / v(s,x)dl s — § / v 2 (s,x)ds ) dx. 



(8.7) 
(8.8) 

(8.9) 
(8.10) 

(8.11) 
(8.12) 
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Since {I t } is an ({J-^}, Q*)-Brownian motion, the conditional expectation simplifies to the 
calculation of a Gaussian integral provided the zero of the max function can be computed. 

With these formulae in hand, we observe that a closed-form expression for the conditional 
expectation can be worked out in the case of a binary initial density function of the form 

fo{x) = qi S(x - x x ) + q 2 6(x- x 2 ). (8.13) 

Here S(x) is the Dirac distribution and q^ = Q[X = Xj\ for % = 1,2. It follows that 



N t = ^Xi q { £ t (xi), A t = ^2 qi £ t (xi), 



(8.14) 



i=l 



8=1 



where we introduce the process 



S t (xi) = exp 



v(s,Xi)dI 8 — ~ I v 2 (s,Xi)ds 



(8.15) 



In the case where the random variable X takes the values Xx and x 2 , the option price is 



c xt = a; 1 e*5 



(8.16) 



We observe that {S(xi)} is positive and has the property E 1 ^ [£ t (xj)] = 1. In particular, we 
can use {S t (xx}} to define a change of measure from Q* to a new measure Q by setting 



(8.17) 



together with dlt = dlt + v(t, x\) dt. Next we pull £t{x\) out to the front of the max function 
in equation (I8.16P to obtain 



c st = A: E 1 



£tM qx{xx - K) + q 2 (x 2 - K) 



£t(x 2 ] 



Lit 



By use of the Bayes formula we express the option price in terms of Q : 

AM 



C st = A; 1 S s (x 1 )E i - 



qx(xx - K) + q 2 (x 2 - K) 



5.19) 



For the sake of a simplified notation we define R ot = S t (x 2 )/S t (xx), and we observe that the 
process {Rot} is an exponential ({J^}, Q)-martingale: 

R ot = exp (^j [v(s,x 2 ) - v(s,xx)]dl s - \ j [v(s, x 2 ) - v(s, xi)] 2 ds^j . (8.20) 

We write Rot = Ro s R s t so that 

C st = A; 1 £ s (xx)E® [( qi (xx -K) + q 2 (x 2 - K)R 0s R st ) + \ F a ] , (8.21) 
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and we note that Rq s is J-^-measurable. One is thus left with the task of finding the range of 
values of R s t for which the max function vanishes. Then we calculate the Gaussian integral 
arising from the conditional expectation. Recalling that {It} is an {{Ft}, Q)-Brownian mo- 
tion, we note that the logarithm of R s t is Gaussian. Let Y be a standard Gaussian variable. 
Then we can write 



ln(i? 



st ) 



[v(u, x 2 ) — v(u, x\)} 2 du Y 



[v(u, x 2 ) — v(u, xi)} 2 du. 



(8.22) 



By solving for the logarithm of R s t in the argument of the max function in (I8.2ip , we deduce 
that the max function is zero for all values y* that Y may take for which 



In 



y < 



q 1 (K-x 1 ) 
q2(K-x 2 )Ro s 



+ \ fsi v ( u i X2 ) - v(u, xi)] 2 du 



'J* [v(u,x 2 ) - v(u } xi)fdu 

It follows therefore that the option price can be written in the form 

1 



(8.23) 



C st = A; 1 £ s (x 1 ] 



q x {x 1 - K) 



. 2vr 

+q 2 (x 2 - K)R 0s 



exp (-\y 2 ) dy 
1 



exp (~\r] 2 {y)) dy 



where 



v(y) = y 



[v(u, x 2 ) — v(u, Xi)} 2 du . 



(8.24) 



(8.25) 



The two Gaussian integrals can be written in terms of the normal distribution function N(x). 
To highlight the similarity with the Black-Scholes option price formula, we define 



d 



St 



-y 



(i 



si 



d st + 



[v(u, x 2 ) — v(u, xi)] 2 du . 



so that one can write 

C st = A; 1 S s ( Xl ) - K) N(d~) + q 2 (x 2 - K)R 0s N(d+) ] . 



(8.26) 



(8.27) 



We can simplify this expression further by use of (I8.14p . Finally, we conclude that the price 
of the call option is given by the following compact formula: 



a 



st 



(xi - K) 



qi + q2 Ros 



N(d; t ) + (x 2 - K) 



q-2 



os + ( h 



N(d 



(8.28) 
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